Graphical prediction of quantum interference-induced transmission nodes 
in functionalized organic molecules 
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Quantum interference (QI) in molecular transport junctions can lead to dramatic reductions of the electron transmission at certain 
energies. In a recent work [Markussen et ai, Nano Lett. 2010, 10, 4260] we showed how the presence of such transmission nodes 
near the Fermi energy can be predicted solely from the structure of a conjugated molecule when the energies of the atomic p z 
orbitals do not vary too much. Here we relax the assumption of equal on-site energies and generalize the graphical scheme 
to molecules containing different atomic species. We use this diagrammatic scheme together with tight-binding and density 
functional theory calculations to investigate QI in linear molecular chains and aromatic molecules with different side groups. For 
the molecular chains we find a linear relation between the position of the transmission nodes and the side group it orbital energy. 
In contrast, the transmission functions of functionalized aromatic molecules generally display a rather complex nodal structure 
due to the interplay between molecular topology and the energy of the side group orbital. 
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1 Introduction 

Quantum interference (QI) effects in molecular junctions has 
recently been suggested as an enabling tool for the imple- 
mentation of molecular switches^ 3 , logic gates 4 , data storage 
elements-^ and thermoelectric devices^— in molecular elec- 
tronics. These concepts originate from mesoscopics, where 
i electron transport through waveguide devices has been in- 
vestigated already two decades ago^~— . In the context of 
single-molecule devices, QI were found to be responsible 
for the observed reduction of the conductance of a benzene 
contacted in the meta configuration as compared to the para 
and ortho configuration s 13 ' 14 , and these findings were ratio- 
nalized by a variety of different physical pictures, such as 
phase shifts of transmission channels or interfering spatial 
pathways^—. More recently, the interest in QI has widened 
to aromatic molecules of increasing sizei^r— and also to in- 
coherent transport in the Coulomb blockade regim e 23 ' 24 . One 
way to induce QI in molecular junctions is to control the elec- 
tron transmission through chemical/conformational modifica- 
tion of side groups to aromatic molecule s Vi 25 " 27 , but also sim- 
pler cross-conjugated molecular wires exhibit QI— and are 
promising candidates for implementing switching and recti- 
fying behaviour—. 

Within the phase-coherent regime, electron transport 
through a molecular junction is described by the energy de- 
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pendent elastic transmission function, T(E). In order to il- 
lustrate what is meant by QI in a molecular junction, it is 
instructive to consider the structure in T(E) as arising from 
three distinct sources. First, we assume that the density of 
states in the electrodes is constant (wide band approximation) 
and that each molecular orbital contributes with an indepen- 
dent channel for electron transport. Under these simplifying 
assumptions, T(E) is a sum of Lorentzian shaped peaks cen- 
tered at the energy of the molecular orbitals (MOs). Next, we 
relax the wide band approximation. This will introduce addi- 
tional structure in T(E); in particular, the peaks in T(E) will 
be shifted and change shape. Finally, we relax the assumption 
of independent transport channels. The structure introduced 
in T(E) in this last step is referred to as a QI effect. The 
impact of QI on T(E) depends on the relative energies and 
shape/symmetry of the MOs, and is in general difficult to pre- 
dict. The most characteristic signature of QI is the presence 
of transmission nodes, i.e. destructive interference, in the en- 
ergy gap between two adjacent MOs. For obvious reasons 
such transmission nodes are most interesting when they ap- 
pear in the gap between the highest occupied molecular orbital 
(HOMO) and lowest unoccupied molecular orbital (LUMO) 
in which case the QI can suppress the conductance by several 
orders of magnitude. 

In a recent work, we showed that the presence/absence of 
QI induced transmission nodes under certain conditions can be 
derived solely from the topological structure of the molecule 
using a very simple graphical methodic. In particular it was 
demonstrated that the graphical scheme correctly predicts QI 
interference in five out of ten different anthraquinone based 
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structures. The graphical scheme is exact for a nearest neigh- 
bour tight-binding (TB) model with equal on-site energies, 
and therefore is expected to work well for all-carbon conju- 
gated systems. The successful application of the scheme to 
systems containing heteroatoms, such as anthraquinone with 
its oxygen side groups, may seem surprising and calls for fur- 
ther analysis. 

In this paper we relax the assumption of equal on-site ener- 
gies. By doing so we can write down equations for the zero 
points of the transmission function directly from the graphs 
leading to a "generalized graphical scheme". For molecules, 
where the conjugated it system is defined by a linear carbon 
chain and a side group (which might contain heteroatoms), the 
scheme predicts a linear relationship between the energy of the 
QI induced transmission node and the energy of the tt orbital 
of the side group. This behavior is confirmed by DFT calcula- 
tions for molecules with a variety of different side groups. For 
aromatic molecules the situation is more complicated due to 
the interplay between topological and on-site energy effects. 
In general several QI induced minima occur in the transmis- 
sion function; their energetic position is given by the roots of 
polynomials in energy which can be derived from the general- 
ized QI graphs. We explore the dependence of these trans- 
mission nodes on the side group on-site energy for the ten 
anthraquinone structures of Ref.— using both TB and DFT 
calculations. 

The paper is organized as follows. In Section [2] we sum- 
marize the graphical QI scheme from Ref.— and generalize it 
to the case of varying on-site energies (heteroatoms). More 
details on the graphical scheme are given in the Appendix. In 
Sections [3] and |4] we present our results for QI in functional- 
ized molecular chains and aromatic molecules, respectively. 
The summary and conclusions are given in Section|5] 

2 Graphical scheme 

In this section we review our graphical approach to quantum 
interference. Some of the details in the derivation are further 
explained in Appendix A. 

Within a single-particle picture, the transmission probabil- 
ity of an electron impinging on a molecular junction with an 
energy E is given by 



T(E) =Tr[Gr i G t r i? ]( J B) 



(1) 



where G = (El — H rao i — — is the Green func- 

tion matrix of the contacted molecule, I is the identity ma- 
trix, Til/r is the self-energy due to the left/right lead, and 

I'l/r = «(Sl/r — Sjr/^). Let us assume that the Hamilto- 
nian describing the molecule is given in terms of a basis con- 
sisting of localized atomic-like orbitals, <fii, <f>2, ■ ■ ■ , 4>n, and 
that only the two orbitals <px and <pjsi couple to the leads. In 



this case the transmission reduces to 



T(E) =j(E) 2 \G 1N (E)\ 



(2) 



Often the energy dependence on the lead coupling strength, 7, 
can be neglected. It then follows that the transport properties 
are entirely governed by the matrix element Gin(E). The 
latter can be obtained using Cramer's rule 



Gin(E) = 



C\n(E1 — H mo i) 



det(£ - H n 



(3) 



where C\n(E — H mo i) is the (IN) co-factor of (El — H mo i) 
defined as the determinant of the matrix obtained by removing 
the 1st row and Ath column from (El — H mo i — S^, — 
and multiplying it by (— 1) 1+N . Since we assume that only 
orbitals 4>i an d 4>n couple to the leads, the removal of the 
1st row and Ath column completely removes ~Sl,r m the co- 
factor. 

In the following we shall focus on C\n(E1 — H mo i) and 
represent the determinant graphically. We use the following 
notation: A hopping matrix element ty , i ^ j is represented 
by a (red) line connecting site i and site j. For simplicity we 
restrict ourself to nearest neighbor hopping, but the applica- 
tion of the graphical scheme is not limited to this case. An 
on-site element (e^ — E) is represented by a (blue) on-site 
loop. Each on-site loop contributes a factor (—1) (see Ap- 
pendix A). We shall set the on-site energy of the "back-bone" 
carbon atoms to zero, £0 = 0, but let the on site energy of 
any side groups, s sg , be of an arbitrary value. This is the main 
difference between this work and the previous worki^S. Also, 
we will not limit the discussion to transmission at the Fermi 
energy, but consider transmission zeros throughout the energy 
range. 

The generalized graphical scheme is summarized as fol- 
lows: The transmission zeros can be determined from the ze- 
ros of the co-factor C\,n(E — H mo i). The terms in the co- 
factor can be represented graphically by drawing all possible 
diagrams according to the rules: (i) In each diagram the ex- 
ternal sites connected to the electrodes must be connected to 
each other by a continuous path, (ii) All internal atomic sites 
must either have one ingoing and one outgoing path or have 
an on-site loop, (iii) The sign of a diagram is (— l) p where p 
is the total number of on-site loops and closed hopping loops. 
When we add up all diagrams constructed according to the 
above rules we obtain a polynomial in E whose roots repre- 
sent the QI induced transmission nodes. 

3 Simple molecular wires 

We first apply our generalized graphical scheme to the simple 
system shown in Fig. [T]where the tt system consists of a nine- 
atom carbon chain (C9) which is cross-conjugated with a side 



2 



1-1 



Fig. 1 Junction setup for a C9 chain with an oxygen side group 
(top). The generalized diagram is shown in the bottom. 

group. Assuming nearest neighbor hopping there is only one 
possible path from left to right which invariably exhibits an 
on-site loop on the side group. This loop corresponds to a term 
— (e sg — E) in the determinant of the cofactor, with e sg being 
the energy of the side group. The path from the left to right 
contacts consists of ten hoppings and thus contributes with a 
factor t 10 , where t is the nearest neighbor hopping energy. To 
find the possible energies at which the transmission is zero, 
i.e. where CiAr(H mo i — EI) = 0, we have only one diagram 
to consider, which is shown in the bottom part of Fig. [U and 
therefore must solve the trivial equation 

- {e sg - E)t 10 = 0^E = e sg , (4) 

for the calculation of the energy, where the transmission zero 
in this case is equal to the on-site energy of the side group. 
We note that the energy of the transmission minimum is in- 
dependent of the on-site energies of the carbon atoms in the 
chain and all hopping parameters in the molecule. Such a 
case where the energy dependence of the transmission min- 
imum is linear with the side group energy will in the follow- 
ing be denoted as a simple side group transmission node. We 
note that such transmission nodes directly caused by the side 
group bears much resemblance with the general concept of a 
Fano resonanc e 31 ' 32 originating from the coupling of a local- 
ized state (the side group) with a continuum (extended states 
in the carbon chain). In Ref.— we have analyzed the analogy 
of the simple side-group transmission nodes with Fano reso- 
nances in more detail. 

Figure |2] (left) shows the transmission for C9-type 
molecules as obtained from DFT calculations^with four dif- 
ferent side groups: CH2, O, NO2, and CF2. For all four 
molecules, we observe a distinct transmission minimum, 
where its energy can be tuned from E — Ep = — 1.4eV for 
O to E - E F = 1.7 eV for N0 2 . In order to illustrate the 
validity of the graphical scheme and the prediction from Eq. 
©, we plot in the right panel of Fig. [2] the energy of the side 
group orbital e sg (see Appendix B for details) vs. the energy 
of the transmission minima Eq. There is a clear linear depen- 
dence (with a slope of 1 .2) in good qualitative agreement with 
Eq. ©. The deviation from a slope of 1 are probably due to 



Fig. 2 Left: DFT calculated transmission functions for C9 with four 
different side groups. Right: tt state eigenenergy of the side groups, 
e 3g , plotted against the energy of the transmission minimum Eq. 
The side group eigenenergy scales linearly with Eq with a slope of 
~1.2. 



the simplifying assumptions in the TB model of only nearest 
neighbour hopping, only a single p z orbital on each atom, etc. 
In the DFT calculations there are interactions beyond nearest 
neighbours and several orbitals on each atom. In addition to 
this more subtle charge transfer effects may lead to deviations 
from the simple TB model. However, overall Eq. (|4|i gives a 
good description of the DFT data. 

4 Aromatic molecules 

We now consider the ten anthracene based molecules illus- 
trated in Fig. |3] which only differ in the position of the two 
side groups. In Ref.— we studied the same molecules with 
oxygen side groups, then they become anthraquinones, and 
categorized them into two groups, where five of them exhib- 
ited QI in the HOMO-LUMO gap (11-15) and the others (Nl- 
N5) did not. We made this distinction on the basis of our 
original graphical scheme and our conclusions were in good 
agreement with DFT calculations. This prediction, however, 
was based on the assumption that the on-site energy of the 
side group would be approximately equal to the on-site en- 
ergy of the carbon p z orbitals (erj ~ Ep). From the side group 
analysis above, we see that this assumption is questionable for 
oxygen with a side group energy of eo ~ — 2eV (relative to 
the Fermi level). We note that this 2 eV difference in on-site 
energies should be considered relative to the size of the hop- 
ping matrix element which is ~ — 3 eV. Still, the success we 
had in applying the graphical scheme^i to the anthraquinone 
molecules deserves a more careful analysis and explanation. 

For this purpose we make use of a single-orbital nearest 
neighbour TB model with carbon on-site energies £0 = and 
hopping elements t = — 2.9eV. We do not explicitly model 
the sulfur end group in the TB model, but include it as an ef- 
fective part of the electrodes. Within this model we vary only 
the on-site energy of the side group and compute the transmis- 
sion function within a wide band approximation, where the 
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Fig. 3 Transmission node energies (black dots) for 10 different anthraquinone molecules as function of the side group on-site energy, e ag , as 
calculated in a nearest neighbour tight-binding model. The dashed curves show the HOMO (lower curve) and LUMO (upper curve) energies. 
The red circles are results from DFT calculations with oxygen side groups. Each panel in the two central columns corresponds to the 
molecular structure in the left- and right most columns. The molecules 11-15 on the left side all exhibit QI in terms of transmission nodes 
within the HOMO-LUMO gap for most values of e 3g , where no nodes are found in this energy range for any of the molecules N1-N5 on the 
right side for any realistic value of e sg - 



results for the transmission zeros are shown in Fig. [3]as black 
dots. The lower and upper dashed lines show the energy of 
the HOMO and LUMO orbitals, respectively, while the red 
circles represent the results of DFT calculations with oxygen 
side groups-^. First, we note that the simple TB model is in 
excellent agreement with the DFT calculations, where both 
result in the same number and the same energetic positions of 



transmission zeros. Our second observation is that the com- 
plexity of the dependence of the transmission minima on the 
on-site energy is quite striking, in particular when compared 
to the simple linear scaling found in the last section for the 
C9-type molecules. 

In order to understand this complexity, we make use of the 
generalized graphical scheme introduced above. However, be- 
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fore turning to the anthraquinones we consider the related but 
simpler benzoquinone structure shown in the upper panel of 
Fig. 2] The lower panel shows the generalized diagrams. In 
addition to the displayed diagrams there are three additional 
but equivalent diagrams where the external sites are connected 
via the lower part of the molecule. These diagrams will sim- 
ply contribute a factor of two to the resulting polynomials. As 
for C9 the on-site loops on the side groups again correspond 
to a term — (e sq — E) in the determinant while for the carbon 
atoms we define — (eo — E) = E with sq = 0. 




Fig. 4 Junction setup for a benzoquinone (top) and all the 
generalized diagrams (bottom). 

Converting the diagrams into a polynomial and taking out 
their common factor (e sg — E)t 3 , we arrive at the equation 

(e sg - E)t 3 [E 2 (e sg - E) - t 2 (e sg - E) + t 2 E] = 0, (5) 

which determines the transmission zeros. We note that it fol- 
lows directly from the diagrams that there will be a maximum 
of four transmission nodes since there can be a maximum of 
four on-site loops no matter how the diagrams are drawn. The 
common factor in the equation corresponds to the common 
part of all three diagrams, where the upper oxygen on-site loop 
contributes (e sg — E) and the path from left to right yields 
t 3 . We thus immediately see that for benzoquinone one of the 
transmission zeros is always at energy E = e sg , which is triv- 
ially defined by the side group on-site energy as for the C9 
molecule in Fig. Q] 

In addition to this simple side group transmission node, Eq. 
§5§ has additional zeros defined by the third order equation 

E 3 - e sg E 2 - 2t 2 E + t 2 e sg = 0. (6) 

While Eq. (|6]l can still be solved analytically, this results in 
rather lengthy expressions for the energies where T = as 



functions of e sg , which do not lend themselves to a transpar- 
ent physical interpretation. Instead we solve the problem nu- 
merically and plot the results for all four transmission zeros in 
Fig. [5] The linear curve (black) corresponds to the simple side 
group transmission node at E — e sg . The remaining three 
zeros exhibit a more complicated behavior, where particularly 
the red curve is of interest because it is situated within the 
HOMO-LUMO gap (marked by the two dashed curves) for the 
whole range of e sg shown in Fig. [5] It is striking that although 
the benzoquinone is much simpler than the anthraquinone 12 
(in terms of the total number of atoms), the nodal structures of 
the two corresponding transmission functions are very similar. 
Most importantly, this applies for both, the simple side group 
transmission node (black curve in Fig. [5]l and the mimimum 
which is situated within the HOMO-LUMO gap (red curve in 
Fig. [5]). Also the remaining two transmission nodes below the 
HOMO and above the LUMO are similar for 12 and benzo- 
quinone. 




Fig. 5 Energy of transmission nodes E, for the benzoquinone as 
function of side group on-site energy, e sg . The dashed lines mark 
the HOMO (lower line) and LUMO (upper line) of the molecule. 
While the simple side group transmission node (black) falls outside 
the HOMO-LUMO gap for large values of |e« 9 |, another node (red 
curve) stays in the HOMO-LUMO gap for all values of e 3g . 

From the five molecules with transmission zeros inside the 
HOMO-LUMO gap, 11-15, only II and 12 have a simple side 
group transmission node. This is straightforward to explain by 
using our generalized graphical scheme. II and 12 are charac- 
terized by having both side groups on the same benzene ring, 
and therefore every path from the left to the right contact has 
to go through at least one carbon atom with a side group di- 
rectly attached. This means that all valid diagrams have at 
least one on-site loop, where (e sg — E) becomes a common 
factor among them, resulting in a simple side group transmis- 
sion node. For 13-15 on the other hand, it is possible to draw 
valid diagrams with no on-site loop on a side group, which we 
illustrate in Fig. |6](top) with a diagram for the central part of 
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13. If such a diagram is possible, no commom factors can be 
found for all diagrams, and hence no simple on-site transmis- 
sion node exists. The nodal spectra of the transmission func- 
tions for 13-15 thus has to be solely defined by higher order 
polynomial equations which we find indeed in Fig. [3] 

We note that all the structures 11-15 have transmission zeros 
within the HOMO-LUMO gap up to rather large side group 
on-site energies |e sg | < 2.0 eV. These are the molecules ex- 
pected to show QI based on the predictions from our original 
graphical schem e 4 ! 30 . This is consistent with the underlying 
assumption of that scheme, namely that the variation in on- 
site energies should be smaller than the hopping matrix ele- 
ment which is ~ 3 eV in the present case. 




Fig. 6 Examples of diagrams for 13 (top) and N3 (bottom). 13 will 
always have on-site loops (blue circles), but not necessarilly on the 
side group. The N3-diagram has no on-site loop and contributes 
with a finite value at all energies. 

We now turn our attention to the five molecules N1-N5, 
which we previously predicted not to have transmission ze- 
ros in the relevant energy range around the Fermi energy 30 . 
Strictly speaking our predictions on the basis of the simpler 
graphical scheme in Ref.— were limited to the assessment of 
the absence or occurence of transmission zeros at the Fermi 
level, when e sg = 0. From the data plotted in Fig. [3] we can 
now confirm that even allowing for a rather large variation of 
the side group on-site energies, namely for |e Sff | < 2.5 eV, 
there are indeed no transmission zeros to be found within the 
HOMO-LUMO gap. We now show how the simple graph- 
ical scheme of Ref.— is a special case of the generalized 
graphical scheme we introduce in the current paper. When 
£sg = £q = 0, all diagrams containing on-site loops will be 
zero, which implies for the 11-15 molecules that all diagrams 
deliver zero terms to the polynomial equations and the trans- 
mission is zero at E = 0. Since there are always two on- 



site loops, as shown in Fig. 6 (top), the transmission node at 
E = eV for e sg = will be doubly degenerate. By chang- 
ing e sg the degenerate nodes split into two and node crossings 
are observed for all 11-15. For the molecules N1-N5 there are, 
however, always diagrams without on-site loops (an example 
is shown for N3 in Fig [6] (bottom)) which will thus contribute 
a finite value to the transmission at E = 0. 

We conclude that the chategorization of the ten 
antraquinone molecules into interfering (11-15) and non- 
interfering (N1-N5) predicted by the simple graphical scheme 
of Ref. 30 holds for all realistic values of e sg . Only for very 
large values of \e sg \, which would require side groups which 
do not couple to the aromatic 7r system and are therefore 
irrelevant for our considerations here, deviations can occur. 
Fig. [3] shows for instance that II has no transmission nodes 
in the HOMO-LUMO gap for |e sg | > 2.5 eV, while for Nl a 
transmission node enters the HOMO-LUMO gap for equally 
rather large |e ag | values. 

5 Conclusions 

We have discussed how characteristic quantum interference 
(QI) induced nodes in the transmission function of conjugated 
molecules can be predicted from simple graphical consider- 
ations involving only the topology of the molecule and the 
on-site energy of non-carbon elements of the 7r-system. A 
previously introduced diagrammatic scheme, strictly valid for 
all-carbon molecules, was shown to be qualitatively valid for 
molecules containing different atomic species as long as the 
on-site energies (p z -orbital energies) do not vary too much 
compared to the interatomic hopping strength - a condition 
we found to be met for a range of conjugated molecules with 
different side groups (O, CH 2 , N0 2 , CF 2 ). We showed that 
more quantitative estimates of the transmission node position 
can be obtained from a straightforward generalization of the 
graphical scheme to the case of finite (and varying) on-site 
energies. This scheme was then used to analyze the transmis- 
sion nodes in linear and aromatic molecules with side groups. 
For linear molecular chains a single transmission node occurs 
at an energy corresponding to the energy of the side group 
7r-orbital while for aromatic molecules the nodal structure of 
the transmission function is in general more complex due to 
a non-trivial interplay between molecular topology and side 
group on site energy. The richness of the nodal structure pro- 
vides a flexible design tool for applications of QI in electronic 
devices based on molecular junctions. 
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A Derivation of the generalized graphical 
scheme 

We shall illustrate the derivation of the graphical scheme by 
considering the four site "molecule" shown in Fig. [7] (top). 
The molecule is connected to the contacts at site 1 and site 4. 
The Hamiltonian describing the molecule is 



tia 

t;i2 




(7) 



El 

tai 


V ui 

As shown in Section [2] transmission zeros are determined by 
the zeros of the co-factor Cin(E — H mo i). The relavant co- 
factor for the four-site molecule is 





£23 
£3 
£43 
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£4 



Ci±{E — H mo i) = 



-*21 E — E2 —t23 

-fa E-e 3 
-ta —£43 



(8) 



The evaluation of the determinant can be done using 
Laplace's formula, 



det(A) 



;(-l) i+j M i: 



(9) 



where the minor Mij, is the determinant of the matrix that 
results from A by removing the i-th row and the j-th column 
(i.e. the cofactor without the sign). Since the on-site ener- 
gies appears in the first upper diagonal at indices + 1), 
each factor (E — £j) have a sign factor (— l) 2l+1 = —1, and 
thus contribute with a minus sign. In the minor l\L l i+ i the 
on-site terms will still be in the first upper diagonal and thus 
contribute with an additional minus sign. 
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Fig. 7 Top: Four site molecule connected to electrodes at site 1 and 
4. Bottom: All the diagrams determining the transmission zeros. 

Writing down the elements of the co-factor gives an equa- 
tion for the transmission zeros (the minus sign is included for 
notational simplicity): 

— C\n{E — H mo i) = 
^12^23^34 - ^14^23^32 + (^2 — E)t 14: (e 3 ~ E) = 0. 



The three terms in the co-factor can be represented graphically 
with the following convention: Each hopping element tij is 
drawn as a line from site i to site j. Note that in all diagrams, 
the terminating sites (in this case site 1 and 4) are connected 
by a continuous path of hopping elements. An on-site term 
(e-i — E) is drawn as a circle around site i and contributes a 
factor (-1). We also note that the loop £23^32 going back and 
forth between site 2 and 3 in the middlemost diagram gives a 
minus sign. It is a general rule, that closed loops, similar to 
on-site loops, contributs a minus sign. This can also be de- 
rived from Laplace's formula: A hopping element tij appears 
at index (i — in the co-factor determinant, Eq. (|8). By 
evaluation of the co-factor along the (i — l)'th row the terms 
involving tij will contribute with a sign (— . The new 
minor appearing after removing the (i — l)'th row and j'th 
column will have the hopping element tji at index (j — 2, i), 
assuming without loss of generality that i > j. In evaluating 
this minor along the j'th column, the term tji comes with a 
sign ( — l) J_2+l and the over all sign of the combination Ujtji 
is (— l) 2l + 2 J -3 = — 1. This shows that a loop between two 
neighbouring sites contributs a minus sign. The over all sign 
of a diagram is thus (— l) p where p is the total number of on- 
site loops and neighbour loops. 

In summary, the transmission zeros can be determined from 
the zeros of the co-factor C\n{E — H mo i). The terms in the 
co-factor can be represented graphically by drawing all pos- 
sible diagrams according to the rules: (i) In each diagram the 
terminal sites connected to the electrodes must be connected 
to each other by a continuous path, (ii) A path can be drawn 
between sites i and j having non-zero hopping elements, tij. 
It is not limited to nearest neighbour interactions only, (iii) 
All remaining, internal sites must either have one ingoing and 
one outgoing path or have an on-site loop, (iv) The sign of a 
diagram is ( — l) p where p is the total number of on-site loops 
and closed hopping loops. 



B Computation of on-site energies 

We calculate the side group on-site energy from the full 
Hamiltonian matrix, H, describing the molecule and the Au 
electrodes. We project onto the subspace spanned by the basis 
functions of the side groups: 



isg 



P HP 

r S griir sg , 



where P sg has diagonal elements on the indices of the side 
group basis functions, and zeros elsewhere. Similarly we get 
a side group overlap matrix, s sg . We then diagonalize s~gh. sg 
to find the side group energies and eigensstates. The corre- 
sponding side group orbitals for the C9 molecules responsible 
for the QI effects are plotted in Fig. [8] 
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Fig. 8 Side group orbitals of C9 responsible for the QI effects close 
to the Fermi energy. 
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